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Abstract. We introduce the discrete Frechet gap and its variants as an 
alternative measure of similarity between polygonal curves. We believe 
that for some applications the new measure (and its variants) may better 
reflect our intuitive notion of similarity than the discrete Frechet distance 
(and its variants), since the latter measure is indifferent to (matched) 
pairs of points that are relatively close to each other. Referring to the 
frogs analogy by which the discrete Frechet distance is often described, 
the discrete Frechet gap is the minimum difference between the longest 
and shortest positions of the leash needed for the frogs to traverse their 
point sequences. 

We present an optimization scheme, which is suitable for any monotone 
function defined for pairs of distances such as the gap and ratio functions. 
We apply this scheme to two variants of the discrete Frechet gap, namely, 
the one-sided discrete Frechet gap with shortcuts and the weak discrete 
Frechet gap, to obtain 0(n 2 log 2 n)-time algorithms for computing them. 


1 Introduction 

Polygonal curves play an important role in many applied areas, such as 3D 
modeling in computer vision, map matching in GIS, and protein backbone 
structural alignment and comparison in computational biology. Given two curves 
in a metric space, it is a challenging task to compare them in a way that will 
reflect our intuitive notion of resemblance. Various similarity measures have been 
suggested and investigated, each of them has its advantages and disadvantages. 
The Frechet distance is a useful and well studied similarity measure that has 
been applied in many fields of research and applications. 

The Frechet distance is often described by an analogy of a man and a dog 
connected by a leash, each walking along a curve from its starting point to its end 
point. Both the man and the dog can control their speed but they are not allowed 
to backtrack. The Frechet distance between the two curves is the minimum length 
of a leash that is sufficient for traversing both curves in this manner. 

Intuitively, the discrete Frechet distance replaces the curves by two sequences 
of points A = (ai,..., a n ) and B = (fq,..., b n ), and replaces the man and dog by 
two frogs (connected by a leash), the A-frog and the i3-frog, initially placed at a\ 
and b\, respectively. At each move, the A-frog or the B -frog (or both) jumps from 
its current point to the next one. The frogs are not allowed to backtrack. We are 



interested in the minimum length of a leash that allows the A-frog and the R-frog 
to reach a n and b n , respectively. The discrete Frechet distance is considered a 
good approximation of the continuous distance, and is easier to compute. 

Much research has been done on the Frechet distance, the majority of which 
considers only the continuous version. However, sometimes the discrete Frechet 
distance is more appropriate. For example, in computational biology where each 
vertex of the polygonal curves represents an alpha-carbon atom. Applying the 
continuous Frechet distance in this case will result in mapping of arbitrary 
points (i.e., interior points on the edges of the curves), which is not meaningful 
biologically. See, e.g., [TJ for a collection of references on the Frechet distance 
and its applications. 

In many of the application domains using the Frechet distance, the curves or 
the sampled sequences of points are generated by physical sensors, such as GPS 
devices. These sensors may generate inaccurate measurements, which we refer to 
as outliers. Several variants of the Frechet distance exist for measuring similarity 
between curves that might be partially erroneous. In particular, variants for 
handling outliers have been proposed, since the Frechet distance is a bottleneck 
(min-max) measure and is very sensitive to outliers. Among these variants are 
the Frechet distance with shortcuts m, the partial Frechet similarity J], and 
the average and summed Frechet distance 03 - 

In the one-sided discrete Frechet distance with shortcuts , we allow the A-frog 
to jump to any point that comes later in its sequence, rather than just to the 
next point. The B frog has to visit all the B points in order, as in the standard 
discrete Frechet distance problem. 

We suggest a new variant of the discrete Frechet distance — the discrete 
Frechet gap. Returning to the frogs analogy, in the discrete Frechet gap the leash 
is elastic and its length is determined by the distance between the frogs. When the 
frogs are at the same location, the length of the leash is zero. The rules governing 
the jumps are the same, i.e., traverse all the points in order, no backtracking. 
We are interested in the minimum gap of the leash, i.e., the minimum difference 
between the longest and shortest positions of the leash needed for the frogs to 
jump from their start points to their end points. 

While the discrete Frechet distance is determined by the (matched) pairs of 
points that are very far from each other and is indifferent towards (matched) 
pairs of points that are very close to each other, the discrete Frechet gap measure 
is sensitive to both. In some cases (though not always), this sensitivity results in 
better reflection of reality; see Figure [l] for examples. 

For handling outliers, we suggest the one-sided discrete Frechet gap with 
shortcuts variant, which we believe has several advantages. Comparing to the 
one-sided discrete Frechet distance with shortcuts, we believe that the gap variant 
better reflects the intuitive notion of resemblance between curves in the presence 
of outliers. Figure [2] depicts two curves that look similar, except for a single 
outlier, with small Frechet gap with shortcuts and large Frechet distance with 
shortcuts. Also notice that the gap variant gives a better matching of the points. 
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(a) (b) (c) 


Fig. 1: Frechet gap vs. Frechet distance: (a) Two non-similar curves, with a large gap. 
(b) Two similar looking curves. The gap is very small while the Frechet distance remains 
the same as in (a), (c) Two non-similar curves, with a small gap and a large Frechet 
distance. 




Fig. 2: (a) The 1-sided Frechet gap with shortcuts is small and the outlier is ignored, 
(b) The 1-sided Frechet distance with shortcuts is large and the outlier is matched. 


Avraham et al. PQ showed that the decision version of the one-sided discrete 
Frechet distance with shortcuts can be solved in linear time, using a greedy 
algorithm. This algorithm can also be used for solving the decision version of 
the gap variant. In this paper, we present an efficient optimization algorithm for 
computing the one-sided discrete Frechet gap with shortcuts, which exploits the 
greediness of the decision algorithm. 

Other variants of the discrete Frechet distance have corresponding meaningful 
gap variants. For example, the weak discrete Frechet distance in which the frogs 
are allowed to jump also backwards to the previous point in their sequence. The 
decision version for this variant can also be solved with a greedy algorithm (in 
quadratic time), and we show how to apply our optimization scheme to efficiently 
compute the weak discrete Frechet gap. In general, our scheme can be applied 
to any variant of the discrete Frechet gap that has an efficient greedy decision 
algorithm. 

Notice that the number of potential gaps is 0(n 4 ), while the number of 
potential distances is only 0{n 2 ). Nevertheless, our algorithms for computing the 
the one-sided discrete Frechet gap with shortcuts and the weak discrete Frechet 
gap run in 0(n 2 log 2 n) time. 

Finally, our scheme can be used for computing the discrete Frechet ratio 
(and its variants), in which we are interested in the minimum ratio between the 
longest and the shortest positions of the leash. More generally, one can replace 
the gap function with any other function g defined for pairs of distances, provided 
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that it is monotone, i.e., for any four distances c < a < b < d, it holds that 
g(a,b) < g(c, d). 


2 Preliminaries 


Let A = (a l5 ..., a n ) and B = (b lt ... ,b n ) be two sequences of points. We define 
a directed graph G = G(A x B,E = Ea U£bU Eab ), whose vertices are all the 
possible positions of the two frogs, and whose edges are all the possible moves 
between positions: E a = j ((a*, bj), {a i+1 , bj)^J j, E B = { ((a*, bj), (a*, b j+1 )^) |, 

Eab = |((aj,6j),(a i+ i,6 J+ i))|. 

The set Ea corresponds to moves where only the A-frog jumps forward, the 
set Eb corresponds to moves where only the .B-frog jumps forward, and the set 
Eab corresponds to moves where both frogs jump forward. Notice that any valid 
sequence of moves of the two frogs (with unlimited leash length) corresponds to 
a path in G from (aq, £»i) to (a„, b n ), and vice versa. 

It is likely that not all positions in A x B are valid; for example, when 
the leash is short. We thus assume that we are given an indicator function 
f:A x B —^ {0,1}, which determines for each position whether it is valid or not. 
Now, we say that a position ( a,i,bj) is a reachable position (w.r.t. /), if there 
exists a path P in G from (ai,6i) to ( at,bj ), consisting of only valid positions, 
i.e., for each position ( ak,bi ) € P, it holds that f(ak,bi) = 1. 

For any distance 6 > 0, the function fg is defined as follows: 


bj ) 


1, bj) < 5 

0, otherwise 


where d(a, : , bj) denotes the Euclidean distance between a, and bj. The discrete 
Frechet distance ddF(A,B) is the smallest S > 0 for which (a n ,b n ) is a reachable 
position w.r.t. f$. 

For any range of distances [s,t], 0 < s < t, the function f[ s t j is defined as 
follows: ( 

1, s < d(ai , bj) < t 
0, otherwise 


f[s,t] ( a i> bj) — 


A range [s,t], s < t, is a feasible (Frechet) range if ( a n ,b n ) is a reachable position 
w.r.t. f [sJ] . 

Let g be a bivariate real function with the following monotonicity property: for 
any four non-negative real numbers c < a < b < d, it holds that g(a , b) < g(c, d). 
Then, the discrete Frechet range w.r.t. g is the smallest 6 = g(s,t) for which 
[s, t] is a feasible range. Two especially important instances are the discrete 
Frechet range w.r.t. t— s, which we name the discrete Frechet gap and denote by 
dFG{A, B), and the discrete Frechet range w.r.t. t/s, which we name the discrete 
Frechet ratio and denote by dpR{A, B). 

In the following sections, we use the difference function t—s as a representative 
function through which we present our ideas and results. However, these results 
are valid for any such function g (and in particular for the function t/s). 
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One-sided shortcuts variants. Let / be an indicator function. We say that a 
position ( ai,bj) is an s-reachable position (w.r.t. /), if there exists a path P in 
G from (ai,&i) to ( ai,bj ), such that /(ai,&i) = 1, /(a,,6j) = 1, and for each 
1 < l < j, there exists a position (ak,h) € P that is valid (i.e., f(ak,bi) = 1). 
We call such a path an s-path. In general, an s-patlr consists of both valid and 
non-valid positions. Consider the path P' (not in G) that is obtained from P by 
deleting the non-valid positions. Then P' corresponds to a sequence of moves of 
the two frogs, where the A-frog is allowed to skip points in A, and with a leash 
satisfying /. Since in any path in G the two indices (of the A-points and of the 
B-points) are monotonically non-decreasing, it follows that in P' the B-frog visits 
each of the points b\,... ,bj, in order, while the A-frog visits only a subset of the 
points aq,..., a,, in order. 

The discrete Frechet distance with one-sided shortcuts d% F (A,B) is the small¬ 
est S > 0 for which ( a ni b n ) is an s-reachable position w.r.t. fs- 

Similarly, a range [s,f] is a feasible (Frechet) range with one-sided shortcuts if 
{a n ,b n ) is an s-reachable position w.r.t. f[ S} t], and the discrete Frechet gap with 
one-sided shortcuts d FG (A, B) is the smallest 5 = t — s > 0 for which [s,f] is a 
feasible range with one-sided shortcuts. 

Weak variants. Let G' = G(A x B , E') be the graph obtained from G by adding 
all backward edges to E, i.e., E' = EU{(v,u)\(u,v) € E}. We say that a position 
( at,bj ) is a w-reachable position (w.r.t. /), if there exists a path P in G' from 
(oi,6i) to ( at,bj ) consisting of only valid positions. Such a path corresponds 
to a sequence of moves of the two frogs, with a leash satisfying /, and when 
backtracking is allowed. 

The weak discrete Frechet distance d^ F (A, B) is the smallest <5 > 0 for which 
(■ a n ,b n ) is a w-reachable position w.r.t. fs- 

Similarly, a range [s,t] is a feasible weak (Frechet) range if ( a n ,b n ) is a w- 
reacliable position w.r.t. /[ s ,t], and the weak discrete Frechet gap d FG (A, B) is 
the smallest 5 = t — s > 0 for which [s, t\ is a feasible weak range. 

3 Computing the discrete Frechet gap 

Observe that for any feasible (Frechet) range [s, t], it holds that t > d dF {A, B) and 
s < min{c?(ai, £»i), d(a n , b n )}. Moreover, since we are interested in the minimum 
feasible range, we may restrict our attention to ranges whose limits are distances 
between points of A and points of B. (Otherwise, we can increase the lower 
limit and decrease the upper limit until they become such ranges.) Thus, we can 
search for the minimum feasible range using the following sorted array of distances: 
D = « in ,..., d f n = min{d(a 1 ,6i), d(a n , &„)},..., d dF (A, B) = ■ ■ ■, C &x ) 

One can compute the discrete Frechet gap d F G{A,B ) in 0(n 4 ) time and 
0(n 2 ) space, by using the standard 0(n 2 ) decision procedure. Start with the 
range [s, t] = [d™ ln , d™ ax ], and as long as the current range is not a feasible range 
increase t (by moving to the next distance in D ). Now, increase s (by moving to 
the next distance in D) and, again, as long as the current range is not a feasible 
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range increase t. Repeat the last step (i.e., increase s by a single move and t by 
a sequence of moves) as much as needed. Finally, return the minimum feasible 
range that was found during this process. 

A recent result of Ben-Avraham, Kaplan and Sharir [2] enables us to reduce 
the running time to 0(n 3 ). Given sequences A and B and an indicator function 
/, they construct a dynamic data structure in 0(n 2 ) time (which also stores the 
information whether (a n ,6 ra ) is a reachable position). Following a single change 
in / (i.e., some valid position becomes non-valid or vice versa), the data structure 
can be updated in 0(n) time. Thus, after increasing s or t by moving to the next 
distance in Z?, we can determine in 0(n) time (instead of 0(n 2 ) time) whether 
[s,t] is a feasible range, since by increasing s a single position becomes non-valid 
and by increasing t a single position becomes valid. 


4 Computing the discrete Frechet gap with one-sided 
shortcuts 

We present an algorithm for computing the discrete Frechet gap with one-sided 
shortcuts in 0(n 2 log 2 n) time. Due to space limitations, the proofs of Lemmas 1-3 
have been moved to Appendix [B| 

4.1 The decision procedure 

Let [s,f], 0 < s < t be a range. We would like to determine whether [s,t] is 
a feasible (Frechet) range with one-sided shortcuts, i.e., whether ( a n ,b n ) is an 
s-reachable position w.r.t. f[ s c ■ In this section, we present a linear-time algorithm 
for doing so, i.e., for solving the decision version. Our algorithm is actually much 
more general and works for any indicator function /. It is similar to the algorithm 
of Ben-Avraham et al. [Tj for the decision version of the discrete Frechet distance 
with one-sided shortcuts. 

Let / be an indicator function, such that /(a*, bj) can be evaluated in constant 
time, for any position ( a,i,bj ). Algorithm [l] computes an s-path in G from (ai,&i) 
to (a n , b n ), if such a path exists. In particular, it determines whether (a„, b n ) is 
an s-reachable position w.r.t. /. Informally, the B-frog jumps forward (one point 
at a time) as long as possible, while the A-frog stays in place, then the A-frog 
makes the smallest forward jump needed to allow the B-frog to continue. The 
frogs continue advancing in this way, until they either reach ( a n , b n ) or get stuck. 

The running time of Algorithm[l]is clearly 0(n ), since the number of iterations 
of the while loop is at most 2n — 1. Notice that we do not construct the graph 
G = G(A x B,E = Ea U Eb U Eab), but the path P produced by the algorithm 
is a path in G. Actually, P is a path in G{A x B,E = Ea U Eb ), since the 
algorithm does not advance the frogs simultaneously. We will use this observation 
later. The correctness of the algorithm is given by the following lemma. 

Lemma 1. Given two sequences of points A = (ai,..., a„) and B = (&i,..., b n ) 
and an indicator function f, dFlS-decision(A,B,f) returns “yes” iff (a n ,b n ) is 
an s-reachable position in G = G{A x B,E = Ea U Eb U Eab) w.r.t. f. 
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Since we did not make any assumptions regarding the function /, Algorithm [l] 
can be used as the decision procedure for both the discrete Frechet distance 
with one-sided shortcuts and the discrete Frechet gap with one-sided shortcuts: 
given a real number 5 > 0 or an interval [s,t], we simply replace / by fs or f[ s t ], 
respectively. 


Algorithm 1 dFlS-decision(A, B , /) 

1. If /(ai,6i) = 0 or f(a n ,b n ) = 0 return “no”. 

2 . curr <— (an,&i). 

3. While true 

Assume curr is ( at,bj ). 

P «— curr. 

If curr is valid (f(curr) = 1) 

If j < n, set curr «— (<n,bj+ 1 ) 

(the B-frog jumps to its next point bj+i, while the A-frog stays at a,i). 
Else ( j = n) return “yes” 

(the A-frog can jump directly to a„, while the B-frog is already at b n ). 
Else ( curr is non-valid, i.e., /(curr) = 0) 

If i < n, set curr <— (oi+i,6j) 

(the A-frog skips a,i, while the B-frog stays at bj). 

Else (i = n) return “no” (the frogs cannot reach position ( a n ,b n )). 


4.2 The search algorithm 


1 2 


Consider the following sorted distances array: 

D = (C in , ■ ■ •, df n = min{d(a l! fer), d(a n , b n )},d s dF (A B) = d r x , ■ • •, C ax ) 
For any feasible range [s,t] with one-sided shortcuts, it holds that t > d dF (A,B) 
and s < min{d(ai, &i), d(a n , b n )}. As in Section[3j we may restrict our attention 
to ranges whose limits are distances between points of A and points of B. 

For the rest of this section, whenever we refer to a feasible range, we actually 
mean a feasible range with one-sided shortcuts. d max 

Let M be the matrix whose rows I-1 

correspond to d 1 ” 111 ,..., d™ ln and whose 
columns correspond to d™ ax ,..., d™ ax (see 
Figure [3]). A cell Mij of the matrix cor¬ 
responds to the range [d™ ln ,d™ ax ]. M is 
sorted in the sense that range con¬ 
tains all the ranges with i' < i,j' < 

j. Thus, we can perform a binary search in 
the middle row to find the smallest feasi¬ 
ble range = [d™ n ,d™ ax ] among the 

ranges in this row. 

Mm. j induces a partition of M into 
4 submatrices: M-\. M- 2 . M :i , M 4 (see Fig¬ 
ure [3]) . Each of the ranges in Mi is con¬ 
tained in a range of the middle row which is not a feasible range, hence none of 



Fig. 3: The matrix of possible gaps. 
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the ranges in Mi is a feasible range. Each of the ranges in M 4 contains 
and hence is at least as large as Mm r Thus, we may ignore M-\ and M 4 and 
focus only on the ranges in the submatrices M 2 and M 3 . 

Sketch of the algorithm. Our goal is to find the smallest range in M for which 
the decision algorithm returns “yes”. This can be done in 0(n 3 log n) time by 
first finding in each of M’s rows (via binary search) the smallest range for which 
the decision algorithm returns “yes”, and then picking the smallest among these 
0(n 2 ) ranges. Below, we sketch a nearly quadratic algorithm for finding the 
smallest feasible range. 

We perform a recursive search in the matrix M. The input to the recursive 
algorithm is a submatrix M' of M and a graph G' by which one can decide for 
each range in M' whether it is a feasible range or not. In each recursive call, we 
perform a binary search in the middle row of M' to find the smallest feasible 
range in this row, using the graph G' . Then, we construct the two graphs for the 
two submatrices of M' in which we still need to search in the next level of the 
recursion. 

Notice that we could use the graph G = G(A x B,E = Ea U Eb) in each of 
the recursive calls, but this would yield an algorithm of running time 0(n 3 logn). 
Instead, in each recursive call we use a graph whose size is proportional to the 
number of rows and columns in the submatrix for this call. The introduction 
of these graphs and their efficient construction is the main contribution of this 
section. 

We represent M and its submatrices by the indices of the array D that 
correspond to the rows and columns of M. For example, we represent M by 
[l,m] x [1, fc]. Mi by [1, f] x [1 ,j - 1], M 2 by [1, f - 1] x [j,k], M 3 by [f + 
1 ,m] x [1 ,j - 1], and M 4 by [f ,m] x \j,k]. 

The skeleton of the algorithm is given below. Recall that g is a bivariate 
real function with the property that for any four non-negative real numbers 
c < a < b < d, it holds that g(a,b) < g{c,d) (see Section [2]). In our case, 
g(s,t) = t- s. 

Algorithm 2 S- Alg(M [Pij/]x [gi9 <], G Mx 

1. Perform a binary search in the middle row of matrix M[ PiP /] x [ 9j9 /] to find the smallest 
feasible range [df™, d™ ax ], i = p , using the decision procedure with the graph 
G[ PlP ']x[g,®']' 

2. Construct the graphs G^i-ijxy.g'] and G[i_|_ ljP '] X [ Sl3 -_i] for the submatrices 

^f[ P ,i_i]xand +1 , p '| x [ry. 1 — 1 ], respectively. 

3. Return 

mm {g{dj ,di ), 

S-Alg(M[ P|i _i] x : G[ Pj i_l] x [j,q'] )i 

S-Alg(M[ i + l ip /] X [q,j — 1] j G [i+iy] x [q,j-l] )}• 


The number of potential feasible ranges is equal to the number of cells in 
M, which is 0(n 4 ). But, since we are looking for the smallest feasible range, 










we do not need to generate all of them. We only use M to illustrate the search 
algorithm, its cells correspond to the potential feasible ranges, but do not contain 
any values. We thus define the size of a submatrix of M by the sum of its number 
of rows and number of columns, for example, M is of size m + k, M 2 is of size 
Y + k — j, and M 3 is of size ^ + j — 1. 

Notice that the ranges in M2 and M 3 consist of all the ranges that intersect 
g = [d™, d“ ax ] and are neither contained in g nor contain g: M 2 consists of all 
the ranges with minimum distance larger than and maximum distance at 
least as large as d™ ax , and M 3 consists of all the ranges with minimum distance 
smaller than df£ n and maximum distance smaller than d™ ax . This implies that 
for any range in M 2 , all the distances in [df l ' lrn , d“ ax ] are in the range and all 
the distances in [ci™, d™] are not in the range (see Figure W|. 
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Fig. 4: The concept of fixed values. 

More generally, let M' be any of the submatrices associated with the i’th level 
of the recursion tree. A distance d £ D is fixed w.r.t. M' if one of the following 
two statements is correct: 


(i) d belongs to all the ranges in M' , or 

(ii) d does not belong to any of the ranges in M'. 


Otherwise, d is non-fixed (w.r.t. M'). The non-fixed distances w.r.t. M' are 
thus the distances that correspond to the rows and to the columns of M' (see 
Figure [4j) . 


Lemma 2. The total size of the matrices in level i of the recursion tree is at 
most m + k, for any level i. 


It follows that the number of non-fixed distances in each level of the recursion 
is at most m + k. We wish to use this fact to reduce the running time of our 
algorithm. 


4.3 The construction of the graphs 

Let dij denote the distance between a* and bj. We say that distance d belongs to 
submatrix M' of M (and write d £ M') if d is one of the distances corresponding 
to the rows and columns of M'. 

We first construct the graph G = G(V = A x B,E = Ea U Eb)- Each vertex 
v i,j = ( ca, bj) £ V has two outgoing edges: 
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1. jumpB(vi t j ) = Vij+i (and if j = n, jumpB(vij) = v n<n ), and 

2 . jumpA(vij) = )_ij (and if j = n, jumpA(vij ) = null). 

Then, we construct the graph Go for the matrix M 0 = M from G, by removing 
the vertices of G whose corresponding distances are fixed w.r.t. M 0 (these are all 
the distances between d™ 111 and d™ ax in the array D), and updating the edge set 
as described below. 

In general, let M,_i and Gi- i, i > 0, be a matrix in level i — l and the graph 
constructed for it; then Vi -1 = {vij \ dij £ Mi- 1 } U {nip, n„ >n }. Let Mi be one 
of the two submatrices of M t _i in level i. We describe how Gi is obtained from 
Gi- 1 (and how Go is obtained from G). Some of the vertices of Gi- 1 are fixed 
w.r.t. Gi, i.e., their corresponding distances are fixed w.r.t. Mi (i.e., they do not 
belong to M,). We say that such a fixed vertex v (whose corresponding distance 
is d) is valid (resp., non-valid) , if d belongs to all ranges in Mi (resp., if d does 
not belong to any of the ranges in Mi). Since backtracking is forbidden, Gi- 1 is 
acyclic and one can topologically sort its vertices. We do so, and then process the 
vertices, one by one, in reverse order (i.e., from last to first). More precisely, for 
each vertex v £ V,_i, we run the code fragment below, where d is the distance 
corresponding to v. The code fragment sets the pointer next[v), for each fixed 
vertex v £ V)_i, so that next(v) is the first non-fixed vertex in the path (in Gi- 1 ) 
beginning at v that is induced by the greedy decision algorithm. 

set next{v ntn ) <- v n ^ n 
if v is fixed and v ^ v n ^ n 

if v is valid (/[ s p](n) = 1 for all ranges [s,t] £ Mi) 
if jumpB(v) is fixed, next(v) 4— next(jumpB(v)) 
if jumpB(v) is not fixed, next(v) ■£- jumpB(y) 

if v is non-valid (/[ s ,t](n) = 0 for all ranges [s, t] £ Mi) 
if jumpA{y) is fixed, next{v) ■£- next(jumpA)v)) 
if jumpA(y) is non-fixed, next(v) £- jumpA(v) 
if v is non-fixed or v = Vi t x 

if jumpB{v) is fixed, jumpB(v) ■£- next(jumpB(v)) (else, do nothing) 

if jumpA{v) is fixed, jumpA{v) ■£- next{jumpA(v)) (else, do nothing) 

Notice that after processing all the vertices of Gi- 1 , it holds that (i) for any 
fixed vertex v v ntn , next{v) is non-fixed, and (ii) for any non-fixed vertex v, 
jumpB(v) and jumpA(v) are also non-fixed (unless maybe when jumpB(v) = 
v n ,n or jumpA{v) = v ntn ). We thus set V) = V)_i\{t> £ V)_i | v is fixed w.r.t. Mi 
and v ui ( i,w n , n } = {vij \ dij £ Mj} U {vi.i,w ni „} and define Gi as the graph 
induced by V). See Figure [5] for an example. 

Given a range [s, t] corresponding to the matrix Mi, we can apply the decision 
procedure to the graph Gp. If both the distances corresponding to v\ t i and to v n ^ n , 
respectively, are within the range (i.e., f[ a ,t]{v i,i) = 1 and f[ s j\(v n ,n) = 1), then 
perform the following loop (otherwise, return “no”). Let v be the current vertex 
(where initially v = nip). If v is valid (i.e., if f[ s .t\(v) = 1), go to jumpB(v ), else 
go to jumpA(v). Return “yes” if and only if you have reached v n , n . This takes 
only 0(| Vi\) time. 
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(b) 


Fig. 5: (a) The graph Go- (b) A graph that was constructed from Go in level 1 of the 
recursion tree. The points represent the vertices of the graphs (all non-fixed). The l’s 
(resp., the 0’s) are the vertices of Go that are fixed valid (resp., fixed non-valid) w.r.t. 
to the new graph. The green arrows mark nextQ, the red ones mark jumpBQ, and the 
blue ones mark jumpA{). The new graph consists only of the points and the red and 
blue arrows. 


Correctness. It remains to prove that the decision obtained when applying the 
decision procedure to Gi is the same as the one obtained when applying it to the 
original graph G. 

Lemma 3. Given a range [s,f] corresponding to M t , the decision algorithm 
applied to Gi returns “yes” if and only if the decision algorithm applied to G 
returns “yes”. 

Running time. Consider the recursion tree. It consists of O(logn) levels, where 
the i’th level is associated with 2* disjoint submatrices of the matrix M. Level 0 
is associated with the matrix Mq = M, level 1 is associated with the submatrices 
M 2 and M 3 of M (see Figure [3|, etc. 

A range test is a test that determines for two distances d± < d 2 in the sorted 
array of distances D whether the range [di,d 2 ] is a feasible range. In the i’th 
level we perform O(logn) range tests in each of the 2 l submatrices associated 
with this level. We claim that the total time spent on the i’th level is 0(n 2 logn). 
This bound includes the preparations towards the next level. Therefore, the 
running time of the entire algorithm is 0(n 2 log 2 n). 

We now focus on the analysis of the i’th level. Let M' be any of the submatrices 
associated with the i’th level. Our algorithm guarantees that the cost of a range 
test, for a range corresponding to M'. is linear in the size of M'. By Lemma [2] 
the total size of the submatrices in level i is m + k < n 2 , and therefore the 
total cost of all range tests performed in the i’th level is 0(n 2 logn). Finally, the 
preparations towards the next level require only 0(?i 2 ) time. 

The following theorem summarizes the main result of this section. 

Theorem 1. Let A = [a\,... ,a n ) and B = (b±,... ,b n ) be two sequences of 
points. Then, the discrete Frechet gap with one-sided shortcuts dp G (A, B) and 
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the discrete Frechet ratio with one-sided shortcuts dp R (A, B) can be computed in 
0(n 2 log 2 n) time. 

5 Computing the weak discrete Frechet gap 

We apply the high-level search algorithm to the weak discrete Frechet gap variant. 
For this we need to (i) describe a suitable greedy decision algorithm and (ii) show 
how to efficiently construct the graphs for the two submatrices of the next level. 
Due to space limitations, we only briefly discuss (i); full details of both (i) and 
(ii) (which is the key issue here) are given in Appendix [A] 

Let G' = G(A x B , E'), where E' = Ea U Eb U {(v, u)|(m, v) G Ea U Eb}- 
That is, G' is obtained from the graph G of the ‘strong’ version, by adding the 
backward edges. (For simplicity, we assume in this section that the frogs are not 
allowed to jump simultaneously.) Notice that G' is a planar graph. We view G' 
as a maze. Each vertex is a room with four doors, one for each outgoing edge, 
that lead to the adjacent rooms (see Figure |6j) . A man standing in room (ai, &i) 
wants to reach room ( a n ,b n ), but without entering forbidden rooms (i.e., rooms 
corresponding to non-valid positions). 

A well known algorithm for traversing a 
maze is the wall-follower rule (also known as 
the right-hand rule): keep your right hand in 
contact with a wall of the maze throughout 
the motion. The algorithm guarantees that you 
will eventually reach the exit, if possible. Thus, 
one can find in 0(nr) time a weak-path in G' 
from (ai,&i) to (a n ,b n ), if such a path exists 
(or reach (ai,&i) if no such path exists). 

The following theorem summarizes the 
main result of this section. 

Theorem 2. Let A = [a\,... ,a n ) and B = (b±,... ,b n ) be two sequences of 
points. Then, the weak discrete Frechet gap dp G (A, B) and the weak discrete 
Frechet ratio dp R (A,B) can be computed in 0(n 2 log 2 n) time. 
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A Computing the weak discrete Frechet gap 


In this section we apply the high-level search algorithm (i.e., Algorithm [2]) to 
the problem of computing the weak discrete Frechet gap. For this we need to 
(i) describe a suitable greedy decision algorithm and (ii) show how to efficiently 
construct the graphs for the two submatrices of the next level. For simplicity, 
we assume in this section that the frogs are not allowed to jump simultaneously, 
however, our solution can be easily adapted to the case where simultaneous jumps 
are allowed. 


A.l The decision procedure 


Let G' = G{A x B, E'), where E' = Ea U E b U {(u, u)\(u, v ) £ Ea U E b }. That 
is, G' is obtained from the graph G of the ‘strong’ version, which contains only 
the forward edges, by adding the backward edges. Let / : A x B — ► {0, 1} be 
an indicator function, which determines for each position whether it is valid or 
not. We provide an 0(n 2 )-algorithm for finding a weak-path in G' (if exists), i.e., 
a path in G' from (ai, b±) to (a„, b n ) that consists of forward and possibly also 
backward edges. 

We describe the weak-path finding algorithm through an analogy to maze 
traversal. Notice that G' is a planar graph. We view G' as a maze. Each vertex 
is a room with four doors, one for each outgoing edge, that lead to the adjacent 
rooms. For a room (a,;, bj). the north and south doors lead to rooms (a»+i,6j) 
and ( a,i-i,bj ), respectively, and the east and west doors lead to (<Zj,6j+i) and 
(a*,6j_i), respectively. All doors are closed, but some are locked and some are 
unlocked. More precisely, a door is unlocked if and only if the rooms on both its 
sides correspond to valid positions. 


—I-1-1-1-1-1— 

0 0 0 1 1 1 1 

_ + + + + + _ l _L r 

1 1 0 1 0 1 n 1 

- + + + + + 


A well known algorithm for travers¬ 
ing a maze is the wall-follower rule 
(also known as the right-hand rule): 
keep your right hand in contact with 
a wall of the maze throughout the mo¬ 
tion. The algorithm guarantees that 
you will eventually reach the exit, if 
possible. (In our setting, whenever our 
hand encounters an unlocked door, we 
go through it as if the wall ended at 
the door.) The wall-follower rule works 
only for simply connected mazes, i.e., 
where all the walls of the maze are con¬ 
nected to the outer boundary of the 

maze, either directly or indirectly through other walls. In our setting, there might 
be rooms that are not connected to the outer boundary (when considering graphs 
for submatrices), but the rule will still work, since the starting and ending points 
are on the boundary. 



Fig. 7: Following the walls (right-hand rule). 
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Thus, starting from the southern wall of room (ai,&i) and using the right- 
hand rule, one can find in 0(n 2 ) time a weak-path in G' from (a\,bi) to ( a n ,b n ), 
if such a path exists (see Figure [?|. If no such path exists, the right-hand rule 
will lead back to (ai,6i), and the decision algorithm will return “no”. 

We note that other algorithms exist for deciding whether a weak-path in G' 
exists (for a given function /), but since this algorithm is greedy, it enables us to 
use the “fixed distances” idea that was used in the previous section. 


A.2 The construction of the graphs 

The notion of fixed and non-fixed distances (introduced in Section [ 4 ]) is relevant 
here as well. Since the rooms in our analogy correspond to vertices in the graph, 
which in turn correspond to distances in D, we simply use the terms: fixed 1-room 
(or fixed 0-room) for a room which corresponds to a distance that is fixed valid (or 
fixed non-valid), and non-fixed room for a room that corresponds to a non-fixed 
distance. 

Initially, we have the graph G'. Each vertex (or room) Vij = (ai,bj) has 
four outgoing (directed) edges (or doors): (i) north(i>ij ) = Vi+ij (and if i = n, 
north(Vij) = null), (ii) ecist(Vij) = Vij + 1 (and if j = n, east(Vij) = null), (iii) 
south(Vij) = Vi-ij (and if i = 1, south(vij) = null), and (iv) west(vij) = 

(and if j = 1, west(vij) = null). 

Now, let Mj_ 1 and Gi_i be a matrix in level i — 1 and the graph constructed 
for it, and let M i: be one of the two submatrices of Mj_i in level i. We describe 
how Gi is obtained from Gi-\. The vertices of Gi correspond to the distances in 
Mi (i.e., we remove from V)_i the vertices that are fixed w.r.t. Mj). 


0 0 



0 0 


1 11111 


+ + + + + G 



Fig. 8: Following the walls (right-hand rule) inside the large room created by a connected 
component of fixed 1-rooms. The missing doors are due to fixed 0-rooms. All the doors 
are to non-fixed rooms. 
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Returning to the maze analogy, we would like to replace the doors leading 
to fixed 0-rooms by walls (we already know they are locked), and to remove the 
doors between pairs of adjacent fixed 1-rooms (we already know they are open). 
By removing doors we get larger rooms with more than four outgoing doors, but 
for each way of entering such a large room, there is exactly one way to exit it, 
using the wall-follower rule. So each incoming edge can be replaced by an edge 
directly to the next non-fixed room. Notice that if the large room has holes in it, 
they can be ignored since the maze’s exit is in the north-east corner and thus 
cannot be in the interior of a large room. 

Gi is constructed from Gi-\ as follows: 

1. Delete all edges (u, v) such that it is a non-fixed room and v is a fixed 0-room. 

2. Let G\_ x be the graph induced by the fixed 1-rooms. For each connected 
component C of G'.p 

— Let Vk.i be any east-most room in C. Clearly, this room is adjacent to 
the outer boundary of C (i.e., to the outer boundary of the union of the 
rooms in C). 

— Starting from the eastern wall of room Vk,i, use the wall-follower rule to 
traverse the outer boundary of C. For each visited room v, let next(v) 
be the first non-fixed room that would follow v if we were walking in the 
graph Gt- 1 (using the wall-follower rule); see Figure [8j 

3. For each edge (it, v) such that u is a non-fixed room and v is a fixed 1-room, 
replace (u,v) by ( u,next(v )). 

4. Finally, delete all the fixed rooms. 

Running time. Computing G' i _ 1 and its connected components can be done 
in 0(|Gj_i|) time using a DFS algorithm. For each connected component C, 
finding its outer boundary can be done in time linear in the size of C, using the 
wall-follower rule. Finding next(v) for each visited fixed 1-room v can be done by 
walking along the outer boundary in the opposite direction. All the other steps 
require time linear in the size of G*_i. We conclude that Gi can be constructed 
in 0(|Gj_i|) time. 

Correctness. The following lemma is analogous to Lemma [6] 

Lemma 4. Given a range [s,i] corresponding to Mi, let be the path traced 

by the decision algorithm in Gj_i. Let Vk and vi be two vertices in 7Ti_i, such 
that i>k and vi are non-fixed w.r.t. Gi but all the vertices between them are fixed 
w.r.t. Gi. Let u be the successor of Vk in 17j_i. Then, 

— if u = north(vk) in Gi-i, then northfvk ) = vi in Gi. 

— if u = east{vk) in Gi- X , then east(vk ) = vi in Gi. 

— if u = south(vk) in Gj_i, then south(vk) = Vi in Gi. 

— if u = west{vk) in Gi- i, then west{vk ) = Vi in Gi. 

Lemma 5. For a given range [s, i], the decision algorithm applied to Gi returns 
“yes” if and only if the decision algorithm applied to G' returns “yes”. 
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Proof. As in Lemma [3j the proof is by induction on i, the level of the recursion. 
Assume that the lemma is true for Gi-±. Let i7i_i be the path traced by the 
decision algorithm in Gj_i. Consider the path 77^ that is obtained from 77,_i 
by removing all the vertices that are fixed w.r.t. Gi. We claim that the decision 
algorithm applied to Gi follows the path 77^. Indeed, for any two consecutive 
vertices Vk and Vk+i in 77^, let u be the successor of Vk in 77i_i. By Lemma[4j 
if u = north{vk ) in Gi -1 then northfvk) = Vk+i in Gi, and the same holds for 
east(vk), south{vk) and west{vk )■ So 77^ is a path in Gi and moreover, since the 
decision algorithm makes only local decisions depending on which direction it 
has currently chosen, it will follow the same path in Gi. 

Now, let v q be the last vertex in 77i_i. If the decision algorithm on Gi -1 
returns “yes”, then v q = v n _ n and the decision algorithm on Gi returns “yes” as 
well. If, on the other hand, the decision algorithm on Gi-\ returns “no”, then 
v q = 1^1 and the decision algorithm on Gi will also return “no”. 

B Missing proofs 

Lemma 1. Given two sequences of points A = (ai, ..., a n ) and B = (bi, ... ,b n ) 
and an indicator function f, dFlS-decision(A, B , /) returns “yes” iff (a n ,b n ) is 
an s-reachable position in G = G{A x B,E = Ea U Eb U Eab) w.r.t. f. 

Proof. It is easy to see that if the algorithm returns “yes” then P is an s-path 
(in G) from (ai,6i) to ( a n ,b n ) and hence ( a n ,b n ) is an s-reachable position in G. 

Assume now that (a n , b n ) is an s-reachable position in G. (Then, in particular, 
f(ai,bi) = 1 and f(a n ,b n ) = 1.) We now prove that if a position v = () 
is an s-reachable position in G, then there exists a position v' = ( ai',bj ) £ P, 
i! < i, such that f(v') = 1. In particular, since (a n , b n ) is an s-reachable position 
in G, there exists a position (a,/ ,b n ) £ P such that f(ai',b n ) = 1, and when this 
position becomes the current position the algorithm returns “yes”. 

We prove this claim by induction on j. The base case where j = 1 is trivial, 
since (ai,&i) £ P. Let 77 be an s-path from (cii,&i) to v = (ai,bj + 1 ). Let 
u = (ak,bj), k < i, be a position in 77 such that f(u) = 1. u is an s-reachable 
position in G, so by the induction hypothesis there exists a vertex v' = (av ,bj) £ 
P, i' < k, such that f(a,i>,bj) = 1. After adding v' to P, the algorithm sets 
curr ■£- (ai',bj+ 1 ). If f(curr) = 1, then we are done. Else f(curr) = 0 and the 
algorithm increases i! until curr = {ai' + i,bj + 1 ), for some l > 1, and f(curr) = 1. 
Since f(v) = 1, we conclude that i' +1 < i, and the claim follows. 

Lemma 2. The total size of the matrices in level i of the recursion tree is at 
most m + k, for any level i. 

Proof. By induction on the level. The only matrix in level 0 is M , and \M\ = m+k. 

Let M 1 be a matrix in level 7 — 1, and assume the size of M' is p + q (it has p rows 

and q columns). In level i we perform a binary search in the middle row of M' 

to find the smallest feasible range [<7™ ln , <7“ ax ] in this row. It is easy to see that 

2 

the resulting two submatrices are of sizes | + q — j and | + j — 1, respectively, 
which sums to p+q — 1. 
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B.l The construction of the graphs in Section [4]- correctness proof 

We assume below that both the distances corresponding to V\% and to u njn , 
respectively, are within the range, since otherwise the claim is clearly true. 

Lemma 6. Given a range [s,i] corresponding to Mi, let Ili-i be the path traced 
by the decision algorithm in Gj_i. Let v k and vi be two vertices in Ili-\, such 
that v k and vi are non-fixed w.r.t. Gi but all the vertices between them are fixed 
w.r.t. Gi. Then in Gi, jumpB(v k ) = vi, if f[ s ,t](v k ) = 1, and jumpA(v k ) = v t , 
if f[s,t](vk) = 0. 

Proof. First, observe that if vi immediately follows Vk in 77,_i, then the lemma 
is clearly true, so let u \,..., u p be the vertices between Vk and vi. We show by 
induction that for any 1 < j < p, next(uj ) = vi. For j = p, if u p is fixed valid, then 
by the decision algorithm vi 4— jumpB(u p ) and by the code fragment next{u p ) = 
jumpB(u p ) = vi. If, on the other hand, u p is fixed non-valid, then by the decision 
algorithm vi = jumpA{u p ) and by the code fragment next{u p ) 4— jumpA(u p ) = 
Vi. For j < p , if Uj is valid, then by the decision algorithm Uj + 1 = jumpB(uf) and 
by the code next(uj ) 4— next(jumpB(uj)) = next(uj+ 1 ), but by the induction 
hypothesis next(uj+ 1 ) = Vi, so we get next('Uj) = Vi. If, on the other hand, Uj is 
non-valid, then by the decision algorithm Uj + \ = jumpA{uj ) and by the code 
and the induction hypothesis next(uj ) 4— next(jumpA(uj )) = next(uj+{) = vi. 

Now, if f[ s ,t]( v k) = 1, then by the decision algorithm u\ = jumpB(v k ) 
and by the code jumpB{vk) 4— next{jumpB{v k )) = next{u\) = vi, and, if 
f[s,t]( v k) = 0, then by the decision algorithm u± = jumpA{vk) and by the code 
jumpA{vk ) 4— next(jumpA(vk)) = next(u\) = vi. 

Notice that Lemma [6] remains true when Vk = wpi, even if v\p is fixed w.r.t. 
Gi, and when vi = v n ^ n , even if is fixed w.r.t. Gi. 

Lemma 3. Given a range [s, i] corresponding to Mi, the decision algorithm 
applied to Gi returns “yes” if and only if the decision algorithm applied to G 
returns “yes”. 

Proof. By induction on i, the level of the recursion. We omit the proof for Gq 
(i.e., that the decision on Go is the same as the one on G), since it is essentially 
identical to the proof of the general case. We thus assume that the lemma is true 
for Gi -1 (i.e., that the decision on Gi-± is the same as the one on G), and prove 
that it is also true for Gj. 

Let IIi_i be the path traced by the decision algorithm in G,_-|. Consider 
the path 77^ that is obtained from 77,_i by removing all the vertices that are 
fixed w.r.t. Gi (except for ui i i j and v n ^ n , even if they are fixed w.r.t. Gi). We 
claim that the decision algorithm applied to Gi follows the path 77^. Indeed, by 
Lemma [bj for any two consecutive vertices Vk and Vk+i in 77^, if f[ s ,t]( v k) = 1 
then jumpB(v k ) = v k+1 and if f[ s ,t](v k ) = 0 then jumpA(v k ) = v k+1 , so 77^ is a 
path in Gi and moreover it is followed by the decision algorithm in G,. 

Now, let v q be the last vertex in Ili-i. If the decision algorithm on Gi-\ 
returns “y es ’\ then v q = v n>n and the decision algorithm on Gi returns “yes” 
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as well. If, on the other hand, the decision algorithm on Gj_i returns “no”, 
then v q ^ v ntn and it holds that either f[ s ,t]{u q ) = 1 and jumpB(v q ) = null or 
f[s,t]i v q ) = 0 and jumpA{v q ) = null. If v q is non-fixed w.r.t. Gj, then by the code 
fragment jumpB(v q ) and jumpA(v q ) do not change and the decision algorithm 
on Gi returns “no” in this case. If v q is fixed w.r.t. Gi, then by the code fragment 
next(v q ) = null. Let vi denote the last vertex in Gi-\ that is non-fixed w.r.t. 
Gi (i.e., vi is the last vertex in ZL). Arguing as in the proof of Lemma [bj we 
get that (in Gi) if f[ s ,t](vi) = 1 then jumpB{vi) = null and if f[ s ,t]{ v l) = 0 then 
jumpA{vi) = null, thus the decision algorithm on Gi returns “no” also in this 
case. 
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